DROPLET MINIMIZERS OF AN ISOPERIMETRIC PROBLEM WITH 
LONG-RANGE INTERACTIONS 
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Abstract. We give a detailed description of the geometry of single droplet patterns in a non- 
local isopcrimctric problem. In particular we focus on the sharp interface limit of the Ohta- 
Kawasaki free energy for diblock copolymers, regarded as a paradigm for those energies modeling 
physical systems characterized by a competition between short and a long-range interactions. 
Exploiting fine properties of the regularity theory for minimal surfaces, we extend previous 
partial results in different directions and give robust tools for the geometric analysis of more 
complex patterns. 



0. Introduction 

In several physical systems competing short-range attractive and long-range repulsive interac- 
tions often lead to the formation of mesoscopic scale patterns. Roughly speaking, the short-range 
interactions favor phase-separation on a microscopic scale, while the long-range ones frustrate such 
an ordering on the scale of the whole sample. When these systems can be described in terms of 
a free energy, such a phenomenon is usually referred to as an energy-driven pattern formation. 
Examples of energy-driven patterns are ubiquitous in physics: among the others we recall ferro- 
magnetic and polymeric systems, type-I superconductor films and Langmuir layers. Even if these 
systems are driven by different physical laws, they exhibit remarkable similarities in the overall 
geometry of the formed patterns (see [33J and [46| ) . 

Our principal interest is the description of the geometry of patterns. For this reason, we leave 
for further studies the detailed analysis of more realistic systems and we focus here on an energy 
model which encodes only the main features of pattern-formation. More specifically, in what 
follows we are interested in the minimization of the following energy functional: 

F 7 , m (u) := / \Du\+j / G(x,y) (u(x) - m)(u(y) - m) dxdy, 
ii» JnJn 

where u is the order parameter of a two-phases system confined in Q C K n , and 7 and m are 
two nonnegativc numerical parameters. The two terms in the energy mimic attractive short-range 
and repulsive long-range energies between the phases. More precisely the first term is local, it 
favors minimal interface area and drives the system toward a partition into few pure phases while 
the second term, involving a Coulomb-like kernel G, is non-local and favor a fine mixing of the 
phases. A detailed description of the energy is given in § [TJ The competition between these two 
terms is expected to induce the formation of highly regular mesoscopic patterns, whose geometry 
strongly depends on the choice of the parameters 7 and m (e.g. spherical spots, cylinders, gyroids 
and lamellae). 

0.1. The Ohta Kawasaki functional for diblock copolymers. The model we consider arises 
as a simplification of a Ginzburg-Landau functional proposed by Ohta and Kawasaki in their 
pioneering paper |41j as a possible description of a diblock copolymers' (DBC) system. Even 
though it is questionable whether such an energy actually describes DBCs (see Choksi and Ren 
[15] . Muratov [39] and Nicthammer and Oshita [40]), nevertheless it is a first, and mathematically 
non-trivial, attempt to capture some of the main features of these systems. For such a reason it 
deserved over the last twenty years great attention from both the mathematical and the physical 
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community (see e.g. [5] [TOl [Ml ED EH] ) ■ Under several simplification, the Otha-Kawasaki functional 
can be written in the following form: 

£ e A u )= / (e 2 |Vu| 2 + W(u)) dx + a / / G(x,y)(u(x) -m)(u(y) -m)dxdy, (0.1) 

where the order parameter u stems for the volume fraction of one block copolymer and W is 
a standard double-well potential. Here m is the average of u over the whole sample, namely 
m = m, and the kernel G solves in Q, the Neumann problem for 

-AG(x,-)=6 x -±- ^G(x,y)dy = 0. (0.2) 

As in the classical Ginzburg-Landau energy functional, the first contribution to the energy forces a 
phase separation through the competition between the gradient and the non-convex potential. On 
the other hand, depending on the strength of the coupling constant a, the long-range contribution 
favors a uniform distribution of the order parameter. This term has an entropic origin in the case 
of DBC (see [8l[l9l|36l|4lJ|49]), but it can also be considered as the energy due to an electrostatic 
interaction between charged bodies if the order parameter is meant to be a density of charges 

(mm). 

It is well-known from the results by Modica and Mortola [38j [37] that, when e -C 1, the 
Ginzburg-Landau energy can be approximated in the sense of T-convergence by a sharp interface 
energy of the form e \Du\, with u being a function of bounded variations taking the two values 
0, 1 (these values identify the pure phases of the system as the sets {x : u{x) = 0} and {x : u(x) = 
1}) and \Du\ denoting the total variation of the measure Du. Formally, this fact gives the link 
between the Ohta-Kawasaki energy and the functional -F 7im (for 7 = a / 'e). It is worth pointing out 
that there exists no rigorous derivation of F lyTn from £ £M in the sense of T-convergence. Indeed, the 
presence of possible multiple scales, (e.g., the one of the phase separation and that of the pattern 
formation) could force the T-limit to be defined on more complex spaces of Young measures as it 
happens in the one dimensional case addressed by Alberti and Miiller [3]. 

0.2. Single droplet minimizers. Physical experiments (see, e.g., [31] and [46]) suggest that, for 
some regimes of the parameters 7 and m, droplets equilibrium configurations are expected. The 
main open issues in this regards are: (1) the rigorous justification of the observed lattice-type 
patterns (for example, the Abrikosov lattice in two dimensions) and (2) the description of the 
geometry of the droplets. 

In this paper we contribute to the second question. In particular, we investigate a regime of 7 
and m leading to the formation of a single droplet minimizer, as a first step towards the analysis 
of multiple droplets patterns. Roughly speaking, a single droplet minimizer can be described as a 
connected region of one phase surrounded by the other one. For this to happen, the competition 
between the two terms of the energy has to be unbalanced, with the confining term stronger than 
the nonlocal one. Moreover, it is also necessary to assume a contribution to the energy coming 
from the interaction with the boundary of f2. This is the reason why in the functional F lt , n the 
total variation of Du is taken in the whole space. If this were not the case, the optimal resulting 
shape would in fact be an almost half ball located in a point of smallest mean curvature of dCl. 
An analysis of this event, though interesting in its own, is not pursued here. 

In order to identify the correct regime leading to a single droplet minimizer, we show here the 
different contributions to the energy of a single ball. As shown in (|1.16j) , given a ball B Tm (p) C ft 
of radius r m centered at p having average mass m, i.e. m|Sl| = uj n r^ (here |f2| stands for the 
n-dimcnsional volume of O), it holds 

2 nr m + 7 (f rt logr + (tt 2 g Tm (p) - ^) r*,) , if n = 2, 
n u n C" 1 + 7 (|^r r?+ 2 + LJ 2 n g Tm (p) r*») , if n > 3, 



F 1 ,m{XB Tm (p)) ~ 
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where g Tm (p) is uniformly bounded for p in a compact subset of f2 - see § 11.41 Therefore, for the 
isoperimctric term to be stronger than the nonlocal one, the regimes to be considered are 

77^ |logr m | << 1 for n = 2, 

jrfn « 1 for n > 3. 

Note that if 7 — > the conditions above are trivially satisfied. On the other hand, in the most 
interesting case of 7 > C > one is forced to consider the small volume-fraction regime r m << 1, 
which we will always assume. Under these scalings we provide a detailed analysis of the minimizers 
of Fj^mi showing that a single droplet is a minimizer for F 1 _ m . In particular, we prove: 

(a) the asymptotic convergence of the minimizers to round spheres in strong norms, providing 
the rate of convergence; 

(b) the asymptotic optimal centering of the droplet in the domain; 

(c) the expansion of the energy in terms of the radius r rn ; 

(d) the nonexistence of exact spherical droplets in domains il different from a ball; and, on the 
other hand, the uniqueness of the minimizer when is a ball (in this case the minimizers 
is itself a ball centered at the center of il). 

These results are summarized in the following theorem (see next sections for more details on 
the notation). 

Theorem 0.1. Let fl C R n be a bounded open set with C 2 boundary. There exist 5o, r o > 
(depending on tt) such that the following holds. A 

7 rf n I log r m \ <S if n = 2 and 7 rf n < 6 if n > 3. 

Then, every minimizer u m = XE m G C m of F^_ rn satisfies the following properties: 

(i) E m is a convex set and there exists p m € fl and ip m : § n_1 — > M such that dE m = 
{Pm + {fm + ^m(^)) x : x € S"" 1 }, for some p rn G fl and ip m : S n_1 — > M and 

\\<p m \\a <7C +3 ; (0.3) 

(ii) p m is close to the set of harmonic centers Ji of£l, i.e. lim rm ^o dist(p m , %) = 0/ 

(iii) the energy of u m has the following asymptotic expansion: 

, x = / 2 7T r m + Y r t lo 8' r m + 7 (- 1 + ^' min n h) rf n + 0(r 6 m ) , n = 2, 
7 ' ml m> \noj n r"- 1 + & C +2 + 7 < r% min fi h + 0(r^+ 2 ), n > 3, 

where h is the Robin function; 

(iv) E m is an exact ball if and only if the domain Q is itself a ball, i.e. up to translations 
Q = Br for some R > 0, in which case E m = B Tm is the unique minimizer. 



Many of the mathematical difficulties in Theorem 10.11 are due to our choice to work in any 
dimension n (previous results are mostly in dimensions n = 2, 3), and with the standard Coulom- 
bian kernel and the natural Neumann boundary condition. Results analogous to ours have been 
obtained under different simplified assumptions in0[T3[Tl[22|32j[39l|42|44l|45]. More in detail, 
in the remarkable paper by Alberti, Choksi and Otto [2] the authors study the uniform distribu- 
tion of the energy and of the order parameter of the minimizers of Fy tm (see [47j for analogous 
results in the case of the Ohta-Kawasaki functional In [39] Muratov studies the shape of the 

minimizers and the asymptotic expansion of the energy in the case of multiple droplet patterns in 
two dimensions for a slightly different nonlocal interaction and periodic boundary conditions. The 
author consider a regime where the minimizers of f 7i?n are given by the union of nearly spherical 
droplets whose centers minimize a pairwise interaction energy (analogous results are also proved 
for the functional £ e ,a)- In a recent preprint by Kniipfer and Muratov [32], the authors study exact 
spherical solutions to a 2-dimensional nonlocal isoperimetric problem in the whole space, where 
the nonlocal term is a Coulombian interaction. In [13l [14] Choksi and Peletier study the regime 
of finitely many (independent of the parameters) droplets proving a T-convergence expansion of 
F 7 m and £ e , a . In (44] [45] Ren and Wei construct explicit equilibria solutions, see § 10.41 for more 
comments. Finally, for a different problem involving a local perturbation of the isoperimetric term, 
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we mention the result by Figalli and Maggi in |22] where an analogous convergence result to single 
droplet minimizers is proved. 

0.3. The role of regularity theory. Many of the existing results regarding the shape of droplets 
are restricted to the two-dimensional case. This is partially due to the fact that in two dimensions 
the isoperimctric confinement is strong enough to allow the use of non-parametric techniques, while 
in higher dimensions having small perimeter does not even imply boundedness (e.g. consider a 
very thin tube). 

One of the main contributions of this paper is to provide robust arguments to overcome this 
difficulty. In particular, we exploit a combination of two facts in the regularity theory of minimal 
surfaces: the uniform regularity properties of minimizers and the use of the optimal quantita- 
tive isoperimetric inequality. More precisely, we show that the minimizers of Fy jTn are uniform 
A-minimizers of the perimeter (see Definition ^. 3j) . thus allowing the uses of non-parametric tech- 
niques. On the other hand, once in a uniform neighborhood of the ball, the optimal quantitative 
isoperimetric inequality is the natural estimate to get the asymptotic convergence for the single 
droplet. We remark that, to our knowledge, this is the first time when the sharp exponent 2 in the 
quantitative isoperimetric inequality has an essential role, while the only case where the uniform 
regularity property plays a role in the same spirit is a recent paper by Acerbi, Fusco and Morini 
[T] where the authors study local minimizers for the functional F im via second variations. 

0.4. Confined solutions to an elliptic system. Problems similar to the ones we have addressed 
here have been considered by Oshita [42] and by Ren and Wei [44] [45] in dimension n = 2, 3. Their 
starting point is somehow different from ours: specifically, they are aimed to the construction of 
special confined solutions to elliptic systems of the form 



Here Q is a bounded smooth domain with boundary d£l having normal v, E C f2 is an open 
set whose boundary dE is a C 2 -manifold embedded in f2, Hqe is its mean curvature and 7 > 
and m £ (0, 1) are two parameters. Note that such a system is the Euler-Lagrangc equation of 
the functional F-y_ rn whenever \e is a smooth critical point. By exploiting a Lyapunov-Schmidt 
reduction procedure, in these papers the authors prove that, for a special range of the parameters 
m and 7, there exists a stationary point of F 1 ^ m having the shape of a single droplet. More 
precisely, they start from a harmonic center p of £1, and by a perturbative argument they are able 
to find solutions to the equilibrium close to the ball centered at p. As a byproduct of Theorem lO.il 
we obtaine the existence of such solutions in any space dimension as the (local) minimizers of the 
associated functional -F 7jm . 

A couple of remarks on this regard are in order. While the analysis made by Oshita is done 
in the somehow simpler regime 7 — > 0, Ren and Wei consider a different range of parameters 
with respect to the one considered here (in particular, they have a restrictive gap condition). 
Surprisingly enough, thought the techniques are different, we obtain the same rate of convergence 
of a single droplet minimizer to the ball in dimension n = 2 as in |44j (while we improve the rate 
in the case of higher dimensions). 

0.5. Variants and possible generalizations. We point out that the techniques developed in 
this paper may also be applied to several other related models, which have been considered in the 
last years. This is because the arguments described in 10.31 do not rely strongly on the form of the 
isoperimetric term neither on that of the nonlocal one, but rather on energy scaling and regularity 
properties of minimizers. Among the possible generalizations, we mention the extension of the 
results of the paper to: 

(1) multiple droplets patterns (work in progress); 

(2) models presenting different Coulombian-type kernels G (such as the screened kernels con- 
sidered by Muratov in [39]): 

(3) droplets minimizers for the Ohta-Kawasaki functional f|0 . 1 [) : 




(0.5) 
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(4) anisotropic perimeter functionals (for which the asymptotic shape will be the Wulff shape 
for the chosen anisotropy - cp. with [2"2"]); 

(5) nonlocal perimeters, as those considered by Carlen ct al. in . 

The paper is organized as follows. In § [T] we fix the notation used throughout the paper and 
recall some known preliminary results which will be used in the proof of Theorem 10.11 In § [5] 
we prove the Lipschitz continuity of the nonlocal part of the energy, deriving the first regularity 
conclusions such as the almost minimality of the minimizers. Then in the short section §|3]we show 
how this observation leads to the main result of this paper in the simpler case when the natural 
Neumann boundary condition are replaced by periodic boundary conditions. We postpone the 
proofs of the general case to § |U In § [5] we discuss the existence of perfectly spherical solutions, 
showing how the regularity plays a role also in the study of the stability. The final Appendix is 
devoted to the proofs of some estimates on the Green function in (|0.2|) we use through the paper. 

1. Notation and preliminaries 

In what follows fi C R n is a bounded open set with C 2 boundary dft. For given constants 
to G (0, 1) and 7 > 0, we consider the sharp interface limit of the Ohta-Kawasaki functional F 1>m 
which can be written in the following way: 

F 1>m (u) := / \Du\+-f / G(x,y)(u(x)-m)(u(y)-m)dxdy. (1.1) 
JR" JnJn 

In the above expression, the order parameter u belongs to the class C m (f2) (often we will simply 
write C m ) of functions with bounded variation taking values in {0, 1}, whose average in fi is to 
and which are constantly outside tt: 



C m := U e BV{R n , {0, 1}) : j u = m, u\ M n\ Q = j . (1.2) 

As already noticed in the introduction, we stress once again that in the functional Ky m the 
total variation of Du is computed in the whole R n , thus accounting also for possible concentration 
of this measure (interfaces of the physical system) on the boundary of f2. In the second term in 
(|1.1[) . G denotes the Green function of the Laplacian with Neumann boundary conditions on dfl. 
Denoting by v the exterior normal to dfl and by \A\ the n-dimensional Lebesgue measure of a set 
A, G is defined by the following boundary value problem: for every i£!!, 

-AG(x, ■) = S x — -p- in O, 
VG(av)-f = on 0fi, (1.3) 

J a G(x,y)dy = 0. 

In place of the average to, we will often make use of the parameter r m corresponding to the 
radius of a ball whose volume fraction in f2 is to, i.e. 

w„r™ :=m|fi|. 

Moreover, we will often identify u € C rn with the set of finite perimeter E such that u = xe (see 
(5J[57j)- We will write the energy F~ m depending on E in the following way: 

F 1>m (E) :=Per(£)+ 7 NL(£;), 

where Pei(E) = J Rn \D\e\ is the perimeter of E in R" and NL is the nonlocal part of the energy. 
Note that, thanks to J n G{x, y) dy = for every x € VL, we may rewrite the nonlocal term as 



NL(£):= / / G{x,y) X E{x)XE{y)dxdy. (1.4) 
Jq Jq 

Finally, we fix the following convention regarding the constants we use in the formula. Every 
time we will use the letter C for a constant, this is assumed to be positive and depending only on 
the dimension of the space n and the domain f2. When possible, we will use the simbols a < b, 
a > b and a ~ b for a < C b, a > C b and C -1 < a < C b, respectively. When we do need to keep 
track of the constants, we will number them accordingly. 
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1.1. Robin function and harmonic centers. Here we recall some basic facts about the Green 
function G. Let T be the fundamental solution of the Laplacian, i.e. 

u n = 2, 



2lT t 2-n (1.5) 

if n > 3, 



r(t) := 

and define the regular part R of the Green function in (jl.3|) as 

: =G(:r,y)+r(|x-y|). 

We recall that, even if in principle R is not defined in y = x, nevertheless, for every x e fi, 
•) solves the following boundary value problem: 

(AR(x,-) = j fl inn, 
\VR(x,-)-u = Vr(\x--\)-v ondfl. 

This implies that R(x, •) is an analytic function in the whole f2 and, hence, it makes sense to 
consider its extension in y = x: 

h(x) := R{x, x). 

The function h is called the Robin function. As it can be easily seen from (|1.6p h turns out to be 
analytic in O. 

Several estimates on the regular part of the Green function and on the Robin function will play 
an important role in the identification of the concentration points for the minimizers of Ey m . The 
following facts will be used in the proofs: there exists r$ depending only on f2 such that, for all 
r <r 

\R(x,y)\~\T(r)\ V x, y : dist(x, dCl) + dist(y, dSl) ~ r. (1.7) 
Moreover, from (jl.7p we deduce also: 

\G(x,y)\<-T(\x-y\) + l, Vx^ef! (1.8) 

h(x) ~ |r(dist(a;,an))| , yxeQ\Q ro , (1.9) 

where, for every r > 0, we denote by il r the complement in f2 of the closed r-neigborhood of 9f2, 
that is 

fl r := {x 6 n : dist(x,dft) > r}. (1.10) 
These estimates are well-known (see for example [24j in the case of Dirichlet boundary condition). 
For readers' convenience we prove these estimate in the Appendix [X] From the regularity of h 
and (|1.9[) . it follows that h is bounded from below. In particular, since h is analytic and blows up 
on the boundary of f2 (hence, in particular it has no constant directions), it follows that the set of 
minimum points of h is finite: we denote this set by % and call them the harmonic centers of f2. 

1.2. The quantitative isoperimetric inequality. The classical isopcrimetric inequality states 
that the perimeter of any measurable set E is bigger than the perimeter of a ball Be having the 
same volume as E: 

Per(.E) - Vev{B E ) > 0, (1.11) 
with equality only in the case E is itself a ball. Quantitative versions of p. lip , also called 
Bonnes en-type inequalities [43] . have been widely studied (see, for instance, [29, 30 ). The follow- 
ing, called sharp quantitative isoperimetric inequality, has been proved in [17( I23|, I25j . 



Proposition 1.1 (Sharp quantitative isoperimetric inequality). There exists a dimensional con- 
stant C = C{n) > such that for every set E C M n of finite measure, it holds 

PerCB^-P^Ss) 
° a{E) ~ P^Be) ' (L12) 

where a(E) is the Frankel asymmetry of E (see [30 ) 

. f f \EA(x + B E )\ 
a(E) = mf | — , i£l 
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Here, VAW = (V \ W) U (W \ V) is the symmetric difference between V and W . For any given 
E C R™ measurable set of positive and finite measure, we say that B ^ is an optimal ball for E if 
\BT\ = \E\ and 

\bT\ 

The center of an optimal ball will also be referred to as an optimal center. In general the optimal 
ball may not be unique but, as proven in |18[ Lemma 6.4], by an elementary application of the 
Brunn-Minkowsky inequality, whenever E is a strictly convex set the optimal ball is actually 
unique. We finally observe that, denoting by r the radius of Be, (|1.12|) scales in r as follows: 

\EAB° E pt \ 2 < r n+1 (Pcr(E) - Pcr(B° pt )). (1.13) 

1.3. First variations. The first variations of F 7>m have been computed for regular sets by Mu- 
ratov |39j in dimension 2 and 3, and then in any dimension by Choksi and Sternberg |16j . Given a 
critical point E of m and x £ dE a regular point of its boundary, the Euler-Lagrange equation 
of F-y. m at E in B r (x) is given by: 

H 9E + 4:jv = c, (1.14) 

where Hqe denotes the scalar mean curvature of dE (namely, Hqe = divfE, with ve the outer 
normal to dE), c 6 R is a constant coming from a Lagrange multiplier and v is the solution of the 
following boundary value problem: 

—Av = xe — rn in f2, 

Vvv = ondn, (1.15) 

Since \\xe — m\\L°° < L it follows that v £ C l a for every a £ (0, 1). Therefore, from standard 
elliptic estimates for the quasilinear mean curvature operator (see [1S])j (|1.14|) implies that, for 
every critical point E, dE is C 3 a for every a <G (0, 1) in a neighborhood of a regular point. As 
shown in the next section, every minimizer of P^ >m is regular except a singular set of Hausdorff 
dimension at most n — 8 (in particular, it is empty in the physical dimensions n = 2,3). 

1.4. Asymptotic energy of balls. Here we give an asymptotic expansion of the energy of small 
round balls in il. Let fl r be defined as in (|1 . 10|) . By the regularity assumption on d£l, there exists 
ro > such that, for every r < ro and p € fl r , the ball B r (jj) £ C Wrjr n. By a direct computation, 

F 7 , Unr n{B r (p)) = Per(B r (p)) + 7 NL(B r (p)) 

= ncj„r™ _1 +7 / / T(\x — y\) dx dy + ■y / / R(x,y)dxdy 



(1.16) 



B r (p)JB r (p) JB r (p)JB r (p) 

= hnr m +j (f rf n logr + (tt 2 g rm (p) - ^) r*) , if n = 2, 
\n Wn r"- 1 + 7 2|^+ 75r (p)( Wll r™) 2 , if n > 3, 

where g r : £l r — > R is given by 

9rip)-=-f f R(x,y)dxdy. (1.17) 

JB r {p) JB r (p) 

In the following lemma we show that g r converges uniformly to the Robin function h as r — > 0. 

Lemma 1.2. Let fl C R" be a bounded open set with C 2 boundary. Then, there exists ro > 
such that, for all r < ro, 

\\g r -h\\ L ~ { n r) ~r 2 , (1.18) 

and, for every r < ro/2, 

g r (p) ~ |r(dist(p,9fi))| y P £fl r \n 2r . (1.19) 
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Proof. To show (|1.18[) . let r be as in (|1.7[) and note that, since -R|a ro xn ro is analytic, we have 
g r {p) - h{p) = -f f (R(x + p,y+ P )-R(p,p))dxdy 

J B r J B r 

{DR(p,p)(x, y) + (D 2 R{p,p) {x, y), (x, y))) dx dy + o(r 2 ) 

D 2 R{p,p) {x, y), (s, y)) dxdy + o(r 2 ), (1.20) 



B r J B 



r 2 



where in the last equality we used that the linear term integrates to 0. By the linearity of the 
integral and of the scalar product, it follows that 

g r ip) - hip) = ^2 {pXid Xj Rip,p) A XiXj + 2 d Xi d yj R(p,p) A XiVj + d yi d yj R(p,p) A yiy ^j, (1.21) 

id 

where 



- 4 — I) 



A-xiXi — A yiyi — fj, :— -j~ x\ dx and A XiXj — A x 

JB X 

By the simmetry i?(x, y) = Riy, x), we infer from (|1.2ip that 

g^p) - h{p) = Tr {D 2 Rip,p)) r 2 + o(r 2 ) = 2/i ARfap) r 2 + o{r 2 ) = + o(r 2 ), 

thus leading to (JTTTHj). The proof of (fTTT9l) follows from (|TT8|) and $L9\). □ 

2. Regularity of minimizers 

In this section we prove the Lipschitz continuity of the nonlocal term, from which we derive the 
uniform regularity properties of the minimizers of F ltTn in the relevant regimes for the parameters 
7 and m. 

2.1. Lipschitz continuity of NL. Proofs of the Lipschitz continuity of NL already appeared in 
the literature (see, for instance, (TJ [39l |48] ) . For our purposes, a more careful quantitative estimate 
of the Lipschitz constant is necessary. 

Proposition 2.1. For every XE m ,XG m G C m , setting w = T * XG m , it holds 

NLiG rn ) - NL(E m ) < (|MU~0E m AG ra ) + \G m \) \E m AG m \. (2.1) 
Proof. We start from (|1.4|) to get 

NL(G m ) - NL(£ m ) = [ f Gix,y)( XGm ix) X G m {y) - XE m {x) X E m (y)) dx dy 

Gix,y)xc m {x) (xG m (y) - XE m iv)) dx dy+ 



QJQ 



QJQ 



+ G { x iV)XE m {y){xG m {x) - XE m ix)) dx dy 



= 2 / G(x, y) xc m (x) (xc m (y) - XE m (y)) dx dy- 



G(x,y) (xG m (y) - XE m {y))- 

ft Jn 

■ {XG m ix) ~ XE m (x)) dxdy, (2.2) 
where in the last equality we used the symmetry G(x, y) = G(y, x). Since 

/ / Gix,y)( XGm {y)-XE m (y))-{xG m ix)-XE m ix)) dx dy = [ |Vz(x)| 2 dx > 0, 
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where z solves 

J ~Az = XG m - XE m in Q, 



Vz ■ v = on dil, 



we deduce: 



NL(G m )-NL(£ ro )<2 / / |G(x,y)|xG m (^)|xG m (2/)-X£ m (?/)| da: dy (2.3) 
in in 

OJ /■ r 

< / / (-T(\ x -V\) + 1 ) XG m ( x ) \xc m (y) - XE m (y)\dx dy 
2 Jn 

(\G m \ - w(y))(xG m (y) ~ XE m {y)) dy 

i 

< (!|w|| L oo (BmAGm) + |G m |) |E m AG m |, 

that is EH). □ 



< 



A straightforward consequence of Proposition 12. II is that, if E m is a minimizer of F i m . then 

Pcr(£ m ) - Pcr(G m ) < 7 (NL(G m ) - NL(£ m )) 

< 7 (lklli~ (£roAGm) + |G m |) \E m AG m \. (2.4) 

By the radial monotonicity of T, it holds 

l|r*x G JU~<l|r*XBj|^= /*/- los ^) ^ = 2. (2 . 5) 



if n > 3, 



k 2 (n-2) 

for every G m with \G m \ = \B rm \. As a result, for r m sufficiently small, we have 

IHU- + |G m |<||r* XG j|^+c<( %^ _l0grm) * n :'' (2.6) 



2 (n-2) 



if n > 3, 



Here we have used the direct computations: 

r * XbM = l*f, + 1 * (1 ° g i r ^ ^ " if - = 2, (2-7) 

\V (logH-i) if \x\>r, 

{ \ x \ 2 i r 2 -f II 

2 " ! - ' if n>3. (2.8) 

—jr. — \ ,„_2 if m > r, 

n (2— n) \x\ n A 117 

Therefore, as soon as r m is small enough that \b t Op) € C m for some p £ fi, it follows by the 
previous two estimates, with G m = B rm (p), that 

mj^)-m^(p))< H'^J 1 ^-Wl * (2.9) 

[77-^ |£? m AB rm (p)| ifn>3. 

In particular, by the quantitative isopcrimetric inequality (|1.13l) . there exists an optimal isopcri- 
metric ball B^ 1 for E m such that 

\E m hB°H | 2 < C +1 (Per(£ m ) - Pcr(S^)) . (2.10) 

In the case XB opt ^m, gathering together (|2.10j) and (|2.9p . we have that 



1 7 r m " n > 3. 



in 
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2.2. Volume constraint. In order to deduce uniform regularity properties for minimizers, it is 
convenient to get rid of the volume constraint. This has been done in several contexts by the means 
of a penalization argument. To this purpose, let us rescale our sets: set p m for the barycenter of 
E m and 

Note that H m is a minimizer of F lr 3^ m in C m (il m ). The following lemma shows that, if H m is 
sufficiently close to a given H C K n well-contained in H m , the volume constraint can be dropped. 

Lemma 2.2. Let H m C fl m be as above for < m < mo with tuq a given constant. Let 
H C £l m C R" be a set of finite perimeter such that dist(iJ, dfl rn ) > 1 for every m €E (0,mo). 
Then, there exists A > with this property: for every m £ (0, mo), if \HAH m \ < A , then H m 
is a minimizer ofG\ im , 

G\,m(E) := F jr s mtTn (E) +A\\E\ — u n \, 
in the class of all sets E with \EAH\ < 2 A -1 . 

The proof of the lemma follows from a simple adaptation of the computations in |21[ Section 2] 
(see also [TJ Proposition 2.7]). We give here only the necessary modifications. 

Proof. The proof is by contradiction. Assume that there exist A; t — > +oo with this property: there 
exist m/i € (0, mo), minimizers of Fh '■= F^ Tm _ rrlh and Eh minimizers of Gh '■= G\ h _ mh such 
that: 

(a) \H h AH\ < A,; 1 ; 

(b) \E h AH\ < 2 A' 1 ; 

(c) \E h \ < \H h \ = uj n (the case \E h \ > \H h \ is analogous); 

(d) G h (E h ) < F h (H h ). 

Since Eh — > H in L 1 (M n ), as in [TJ Proposition 2.7], one can show the existence of suitable 
deformations Eh satisfying the following: 

(e) |^ h | = \H h \=cj n ; 

(f) dist(E h ,H) < 1 (in particular, E h C O m ); 

(g) there exist 07,. > with \Hh\ — \Eh\ > Ci(n) ah such that 

Per(^) < Per(S h ) (1 + c 2 (n) a h ) and |f h A£? fc | < c 3 (n) a h Per(E h ), 
where ci, 02,03 > are dimensional constants. 
Hence, we infer that, for h sufficiently large, 

F h (E h ) = G h (E h ) 

< G h {E h ) + [c 2 {n) G h Vcy{E h ) + C \E h AE h \ -A h \\E h \- w n |] 
C^).(ff) 

< F h (H h ) + a h [ca Per(E h ) + C c 3 Per( J B 7l ) - ci A h ] 

< F h (H h ), 

where we used the Lipschitz continuity of the nonlocal term, Ah — > +00 and the uniform bound 
on Pei(Eh) implied by: 

Pcr(E h ) < F h {H h ) < F jr 3 , mh ( B i) < +°° V ft G N. 

h 

This contradicts the minimizing property of Hh and, hence, proves the lemma. □ 

2.3. A-minimizers. It follows from Lemma |2~21 that the sets H m are uniform strong A-minimizer 
of the perimeter according to the following definition. 

Definition 2.3. Let fl C K™ be open. A set of finite perimeter E C f2 is a strong A-minimizer of 
the perimeter in f2 at scale r\ > if 

P{E) < P{F) + A\EAF\, V EAF CC n, \EAF\ < 77. (2.12) 
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This notion of almost minimality has been widely studied in the regularity theory for minimal 
surfaces. By the theory developed in [?J |H1 HH HO] (in particular, see [SOI Theorem 1.9 and Propo- 
sition 3.4]), strong A-minimizcrs have regularity estimates which arc uniform in the parameters A 
and r\. More precisely, for every a £ (0, 1), there exists a constant e = £q(ti, a, A, rj) such that 

Exc(£;,B r (x)) := r 1 "" (Pcr(£) - \D X E(B r (x))\) < e a => dE n B*(x) £ C l a . 

Since the quantity Exc is continuous under L 1 convergence of A-minimizers, uniform regularity 
estimates can be inferred for A-minimizers in a neighborhood of a given smooth set. 

Proposition 2.4. Let E d Q be a strong A-minimizers at scale n and let F C R" be a set with 
smooth boundary and dist(i 7 ', dfl) > 1. Then, for every a £ (0,1), there exist constants r)o = 
r)o(n, a, A, rj) > 0, R — R(n,A,r)) > 0, c = c(n) > and a modulus of continuity u> : R + — ¥ R + 
with this property: 

(i) if \EAF\ < TjO) iften 9-B can 6e parametrized on dF by a function if : dF — > R. 

dE = {x + <p(x) v F (x) : x £ dF}, 

with IMIc 1 ^ < u(\EAF\); 

(ii) for all x £ E and < r < R with B r (x) C f2, it holds 

c(n) r" < \EnB r (x)\. (2.13) 
Note that, in the sequel, we will apply Proposition 12.41 always in the case F = B\. 



2.4. Higher regularity. Thanks to Proposition 12. 4[ the first variations in § 11.31 can be used to 
improve the regularity of the minimizers of F 1 ^ m . 



Proposition 2.5. Let E m be a minimizer of F 1<m and let H m , p rn and Q m be as in § \2.2l Then, 
for every a £ (0, 1), there exists r\ > such that, if^yr^ < 1, \H m AB\ | < rj and dist(_Bi, <9f2 m ) > 1, 
then H m can be parametrized on dB±, 

8H m = {(1 + ip m {x))x : x £ dBx}, 

and \\<Pm\\c 3 < a ^ &{\H m ABi\) for a given modulus of continuity us . 

Proof. The existence of a parametrization (p m is guaranteed by Proposition 12.41 (i), under the 
hypothesis that n is chosen sufficiently small. We need only to show that the Euler-Lagrange 
equation for F 1 _ rn , namely 

HdH m {x + ip m (x) x) +4^rf H w m (x + (p m (x)x) = A m , (2-14) 

with A m £ R a Lagrange multiplier, allows actually to infer the higher regularity claimed in the 
statement. Note that it suffices to prove sup TO ^ H^mllc 3 ." < C (note that the parameter 7 is 
hidden in the notation tp m ). Indeed, since HVtoIIc 1 -" ^ ^{\H m AB\\) — ¥ as rj — > 0, where uj is 
the modulus of continuity in Proposition 12. 4| by compactness in the C 3 ' Q norm we would as well 
deduce that ||<p rn || C 3,o — > 0. 

To show this, we consider separately the two different terms in (|2 . 14[) . For what concerns 
X rn we recall that, by Lemma 12.21 there exists > such that H m minimize Ge,m locally in a 
neighborhood of B\. This allows us to compute the first variations of Ge.m- Since the penalization 
term \\E\ — tu n \ is not differentiable, we have to distinguish between the variations increasing 
the volume and those decreasing it. Let ip £ C°°(dBi) and K e be the competitor 

dK e := {x + ((p m (x) + e ip(x)) x : x £ dBi}. 

The volume of K e is given by 



\K e \=u n (l + ip m + etp) n dH n -\ 

JdB x 

hence, it follows that \K e \ > (J n or \K e \ < u> n for small e>0if / i/>>0or J ijj < 0, respectively. 
The minimizing property of H m implies the following variational inequality to hold true: 

Je le = 0+ -°- 
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In turns this leads to (with analogous computations for the first variations of F~ m as in |16j ) 
(H dHm (x + (p m {x)x)+Ayrl l w m (x + ip m {x)x)+e^-ip{x)>0 if / i>dU n - l >0, 



(2.15) 



< 0, 



[HdH m {x + (p m {x)x) + 4jrf n w m (x + ip m (x)x) - 0) tp(x) > if / ipdH r ' 

(2.16) 



where w m solves the boundary value problem: 

-Aw m = XH m - m in Q m , 

Vw m ■ v = on dtt m , (2-17) 

Jn m w ™ = °- 

Hence, from (|2.15p and (|2.16|) we deduce a uniform bound on the Lagrange multipliers A m : 

|A m |<6 y m >0. (2.18) 

For what concerns w m , by an analogous computation as in (|2.5[) using \G\ < |T| + 1 and the radial 
monotonicity of T, we deduce that ||w m |jL°° < C. Moreover, since \\xH m — Xb x \\lp ^ '/ for every 
p > n, the Sobolev embendings and the Gagliardo-Niremberg interpolation inequality leads to 
uniform W 2,p bounds and, hence, C 1,a bounds on w m for every a £ (0, 1) (see Chapter 9]). 

Hence, since (p m has also uniform C > a bounds, the non-parametric theory for the mean 
curvature- type equation (|2.14j) (see [33J Chapter 3] or [571 Appendix C]) finally yelds the desired 
uniform C 3 ' a estimates for <p m . □ 

3. Periodic boundary conditions: fl = T" 

In this section, in order to show the proof of our main result in a technically simpler case, we 
prove a statement analogous to Theorem 10.11 for periodic boundary conditions. Indeed, since in 
the present case one discards the interactions with the boundary and the optimal centering of the 
asymptotic droplet, the proof is a direct consequence of the regularity arguments developed in the 
previous section. 

3.1. Notation and statement. Let T n be the n-dimensional torus obtained as the quotient of 
K n via the Z n lattice or, equivalently, T™ := [0, 1]™ with the identification of opposite faces. We 
consider functions 



£ BV(T n ;{0,l}) with f 

JT" 



As usual such functions u can be identified with measurable sets E C W 1 invariant under the 
action of Z™ and such that \E H [0, 1]"| = m. The confining term of the energy is then given by 
the perimeter of E in the torus: 



Per(£,T") := / \D X 

y[o,i)» 

and the nonlocal term by: 



NL(£):=/ / G{x,y)xE(x) X E(y)dxdy, 

J[0,l] n J[0,l] n 

where G is the Green function for the Laplacian in T™, i.e. 

f-AG(av) =5,-1 inT", 
\j Q G(x,y)dy = 0. 

By the invariance of the torus, we can write with a sligth abuse of notation G(x, y) = G(\x — y\). 
In the case of periodic boundary conditions, Theorem 10.11 reduces to a statement regarding the 
shape of the minimizers E m and the asymptotic behavior of the energy. 
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Theorem 3.1. There exists S > such that the following holds. Assume r m < 1 and 

7 r^|logr m | < S if n = 2 and jr 3 m < S if n > 3. 
Then, every E m C T" minimizer of Ej jTn is, up to a translation, a convex set such that 

dE m = {(1 + i> m (x)) r m x:xe S"" 1 }, 
for some ip m : IS" -1 — > R with 

Mm\\ci<ir%+ 3 , (3.2) 
and its energy has the following asymptotic expansion: 

F ( , = /27rr m + ^logr m + 7 K + 7r 2 M0)) r* + 0(r*J, n = 2, 
" MXEm) k^C- 1 + |^fC+ 2 + 7 ^ r ^M0) + O(e +2 ), n>3. 1 • > 



3.2. Improved perimeter estimate. Due to the translation invariance, for every minimizer E m 
we may assume that B ^ 1 = B Tm is centered at the origin. Therefore, from p. lip we infer for 
H m = E m /r m that 

iff a r I < J 7r ™ I lo S r '"l < 5 ° if n = 2, 
^7 r m < if n > 3. 

If Jo is chosen sufficiently small, by Lemma 12.21 the sets H m are minimizer of some penalized 
functional G^ m and, hence, are A-minimizers of the perimeter. By Proposition 12.41 H rn can be 
parametrized by the graph of a function tp m on dB\ satisfying 

||¥>ml|l,°°(8Bi) ~ \H m /\B±\, 
thus implying that E m can be parametrized on dB Tm by ip m with 

II , II <- l-^mAB^J 

||Vm||i,°°(8B rm ) £ — ■ (3-4) 

These observations lead to the following proposition which is a consequence of an improved esti- 
mate for the Lipschitz constant of the nonlocal part of the energy. 

Proposition 3.2. There exists 5q > such that, if 7 r^j |logr m | < Sq in the case n = 2 and if 
7?"m < So in the case n > 3, and E m is a minimizer of F 1%rn , then 

\E AB opt \ 2 

?er(E m ) - Per(B^) < 7 ' " »_f J + 7 C +1 \E m AB°£\. (3.5) 
Proof. Recalling (|2.3p in Proposition 12.1) and assuming as above B rm = B°g , we have that 
NL(S rm ) - NL(£ m ) < ( / G(x, y) X B Tm (x) ( X B rm (y) - XE m (y)) dx dy 

r(k - y\)xB Tm (x) (xB rm (y) -XE m (y))dxdy+ 
R(x, y) XB rm (x) (xB rm (y) - XE m (y)) dx dy. 

nJn 

By the direct computation of w = T * \B r in l|2.7[) and (|2.8[) (in particular, |Vio| < r m in a 
neighborhood of dB Tm ), we get, 

|£! m AB r I 

HHU°°(S m AB rm ) ^ r rn ||^m||i«>(8B rm ) ^ 



n Jn 



r n-2 
1 m 



Moreover, 



Q JQ 



R{x,y) XB Tm (x) (xB rm (y) - XE m (y)) dx dy 



= #(0,0)/ / XB rm (x)(xB rm (y)-XE m (y))dxdy + OK+ 1 )\E m AB r J 
Jn Jn 



r 



™ +1 \E rn AB ri 
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where we used that J Xb t ~ J XE m = 0- Hence, gathering together the previous estimates, by 
the minimality of E m , it follows: 

Per(i; m ) - Pcv(B rm ) < yNL(B r J - 7 NL(£ m ) 

j\H\L°°(E m AB rm ) \E m AB r J +<yr^ 1 \E m AB rm \ 

^l lEm t B 2 rJ2 +ir n m +1 \E m AB r J. 

□ 

3.3. Proof of Theorem 13.11 In order to prove (|3.2[) we use the quantitative isoperimetric in- 
equality and the improved estimate in Proposition 13 . 21 to infer that 



\E m AB°£\ 2 < C+^Per^) - Per(B^)) 

< ir 3 m \E m ABTj 2 +irl^ \E m AB°£\. 
This implies \E m AB^\ < 7r 2 a ™+ 2 , that is, by flO}, 

!lVAn||L 1( a Bl )<7< +3 - (3-6) 

From the C 3 a regularity of ijj m proved in Proposition 12.51 the convexity of E m and (|3.2[) follows. 
Similarly, by comparing the energy of E m with that of B Tm , using Proposition ^. 1[ Proposition ^. 21 
and Lemma fOl (|3 . 3[) follows: 

F y , m (E m ) = F 7 , m (B rm ) + 0( 7 ^™+ 3 ) 

= hnr m + ri logr m +7 (-| + tt 2 fc(0)) r r 4 n + 0(r£j, if n = 2 

" \noj n r?- 1 + r™+ 2 + 7 c 2 r 2 ? fc(0) + 0(rf™+ 2 ), if n > 3. □ 

4. Strong convergence to round spheres 
In this section we prove Theorem lO.il (i). (ii) and (hi). 

We remark that, in this general case, before we may argue as in the proof of Theorem 13. li we 
need to show that the minimizers of i^ )TO are well-contained in 0. This is crucial in order to apply 
the regularity results in § [2j which hold under the hypothesis of being at a fixed distance to the 
boundary. The proof is based on the analysis of the non local energy of a minimizer when it gets 
close to dfl. To this extent a key role will be played by the estimates (|1.7[) . (|1.8[) and p.9[) . 

4.1. Localization of minimizers. In the next proposition we prove that the minimizers of F 1<m 
are well-contained in Q. 

Proposition 4.1. There exist So,ro > such that the following holds. Assume r m < ro/3 and 
"fr^ |logr m | < Sq ifn = 2 and 7r 3 „ < So ifn > 3. Then, every minimizer E m of F 7)TO satisfies 

E m C B 3rm (q) for some q £ f2 ro . (4.1) 

Proof. We prove the result in the case n > 3, the case n = 2 being similar up to minor changes. 
The proof consists of three steps. 

Step 1. If Sq and ro are sufficiently small, then there exists a ball B m := B rm (p m ) C f2 such 
that 

\B m AE m \ <5q /( " +1) C and dist( Pm ,dn) > S- 1/{{n+1){n ~ 2)) r m . (4.2) 

For any ball of radius B rm (p) C O (note that such a ball exists if ro is chooscn sufficiently 
small), by flUJ it holds 

Pev(E m ) - Pev(B rm (p)) < -/r 2 m \E m AB rm (p)\ < 7C +2 - (4-3) 
On the other hand, by the quantitative isoperimetric inequality we have 



\E rn AB%/ ™ C+ 1 (Per(£ ro ) - Per(fl^)) T 1^ +3 . (4.4) 
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Note that B° P * may not be contained in Q. Nevertheless, since Btf \ n C B^AEm, it follows 
from (|4~4f that 

l^\n|<(7d)Va C <6l'\" m . 

We can now use a simple geometric argument proved in Lemma 14.31 below to deduce the existence 
of a vector v £ M. n such that 

M ~ S 1 /{n+1) r m and B m := + v C O. 

Setting p m as the center of B m , namely B m = B rm (p m ), this amounts to say that p m €E ^r m - We 
claim that B m satisfies (|4.2j) . Note that, since the measure of the symmetric difference between 
two balls is linear with the distance of the centers, we infer the first conclusion in (j4.2[) . namely 

\E m AB m \ < \E m AB°£l\ + \Bf n AB m \ < (<^ /2 + r« < <5 1/( " +1) C. (4.5) 

On the other hand, appealing to the minimality of E m (now XB m € C m ) and using (|1.16p , we get: 
7(cj n r™) 2 g rm {p7n) - 7(^n<J 2 min g Vm (p) = F im {B m ) - min Fy )m (B rm (p)) 

<7iVL(S m )- 7 AT J L(S m ) 
<7^o 1/( " +1) C +2 , (4-6) 



where in the last inequality we have used Proposition 12.11 and (|2.6I) . Then, by Lemma 11.21 and 
3]) we obtain the second inequality in 



{Ai^p^dn) +r m ) 2 ~ n <&y {n+1) r* m n + min g Tm (p) < 

Step 2. TTie whole E m is well- contained in f2, i.e. 

E m d B 3rm (p m ). (4.7) 

With the notation as in S 12. 2[ by (|2.11|) we have that \H m ABi\ < 77^ < 5 . Then, using 
Lemma 12. 2\ the sequence of sets H m turns out to be a sequence of uniform A-minimizer of the 
perimeter in Q m . Moreover, by (|4.2jl . if 5q is small enough, we have that 

dist(o, dn m ) > ^ o -v(»+i)(«-2)) > 4 . (48) 

As a consequence, we are in position to use the density estimate (|2.13|) in Proposition ^. 41 according 
to which there exists R > (without loss of generality we assume R < 1) such that, for every 

xeH m n(B 3 \B 2 ), 

c(n) iT < |# m n < \H m n B!(a:)|. 

Therefore, since for every x G H m D (S3 \ -B2) it holds £?i(x) n Si = 0, we get: 

c(n)R n < |H m nBi(i)| < |ff m ASi| < <5 1/(n+1) . 

Clearly, if 5 is small enough, this inequality cannot be satisfied, thus implying H m n (B 3 \B 2 ) = 0. 
In order to complete the proof of (|4.7[) . we need to show that H m n (f2 m \ B 3 ) = as well. To 
this purpose, we reason by contradiction and show that, in this eventuality, a suitable rescaling of 
J m := H m n B 2 would have lower energy than H m . We fix the notation: 

A ni • — B- m \ Jyyi and Lm ■ — Pm J m , 

with p m > 1 such that \L m \ = \H m \. Note first the following two observations: (a) by a simple 
computation on the volumes, it follows that 

p m -l<\K m \; (4.9) 
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(b) consequently, we can estimate \L m AJ m \ in the following way: 
\L m AJ m \ = I \XJ m (Pm x ) - XJ m ( x )\ dx 



< / / \D X j m (sx + (l-s)p- 1 x)\(l- p^)\x\dsdx 
Jb 3 Jo 

<(p m -l)Per(J m )<\K m \, (4.10) 

where, in order to rigourously justify the second inequality without referring to fine properties of 
functions of bounded variations, it is enough to consider an approximation with smooth functions 
and to pass to the limit. Recalling that F^ tm (E m ) = r™~ 1 i ? 7r 3 m (H m ), we can compare the 
energies of H m and J m as follows: 

F jr 3 m ^ n (L m ) = p^ 1 Per(J m ) +7^ NL{L m ) 

< p^- 1 Per(F m ) - p™- 1 Per(A m ) + 7 r m C \L m AH rn \ + 7 r 3 m NL(H m ) 

< 9 (1 + C \K m \)Per(H m ) - p^ 1 Pex(K m )+ 
+ 7 r 3 m C (\L m A J m \ + \K m \) + 7 r 3 m NL(H m ) 

^F ir 3 m . m (H m ) + C \K m \ + CS Q \K m \ - C |A m | ( "- 1)/n 

< F jr s m!Un (H m ), (4.11) 

if Sq is sufficiently small since 

I if ml < Sl /{n+1) < 1. Clearly, this is a contradiction with the 
minimality of H m , thus proving that H m \ B3 = or, after scaling by r m , that (|4.7|) holds true. 

Step 3. Proof of (|4.1[) . We set E' m := E m — p m and, as a consequence of (|4.7[) . wc note 
that E' m C B 3rm . For all q £ ^3r m , let us set E m (q) := E' m + q (in particular, E m {q) C f2 and 
E m = E m {p m )). We may write the energy of E m (q) as 



F 7im (E(q))=Pei(E'J+ 1 / / T(\x - y\) X E> m (x) X E> m (y) dx dy + 

+ 7 / / R(x + q,y + q)xE'(x)XE'(y)dxdy. (4.12) 



Since E m minimizes -F 7 ,m, we have that F 7 iTn (E m (p m )) < F 7 ^ m (E m (q)) for every q £ ^3r m - By 
(|4.12j) this implies that 



// 



R(x - p m , y - p m ) XE' m (x) XE' m (y) dx dy 

< / / R(x-q,y-q)xE'(x)xE'(y)dxdy. (4.13) 



In view of E' m C £?3r m , (|1 -7p and (|1.5[) the last inequaltiy (|4. 13[) implies that p m is contained in a 
compact subset of fi, thus proving (|4.1[) . □ 

Remark 4.2. It follows a posteriori that the optimal balls B ^ for E m are, in fact, well-contained 
in D, and (14.41) holds, i.e. 



dist(B^,afi) > a-VKn+Wn-a)) and \B^AE m \ < <5o /2 r™, (4.14) 

The following is the geometric lemma used in the localization argument. 

Lemma 4.3. Let f2 C K n be an open set with C 2 boundary. Then, there exist ro, ho > with this 
property: for r < ro, h < ho and p £ f2 such that \B r (p) \ tt\ < hr n , there exists v £ R™ with 
\v\ < /i 2 /(" +1 )r such that B r (p + v)cn. 
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Proof. The main argument in the proof is given by an elementary consideration. Assume first 
that dfl n B x (p) is flat. If \Bi(p) \ 0| < h and h < h is small enough, then dist(g,r) ~ h 2 /( n+1 \ 
To see this, one can easily compute the exact expression for (3 := dist(g r , r) solving the equation 



h = (n-l)u n -ij [y/l-r 2 - 1 + /3J r n ~ 2 da;. 

Alternatively, one can simply notice that dist(g, s) ~ /3 1 / 2 and the volume of B\ (p)\f2 is comparable 
with that of the cylinder with base dil n B 1 (p) and height /3 (in fact, the cylinder with half the 
height and half the radius of the base is contained in B\{p) \ ft). Hence, h ~ f3( n + 1 ')/ 2 from which 
the conclusion. Clearly, v — —(3 e n fulfills the conclusion of the lemma. 

If f2 is not flat, we need to restrict the size of the balls we consider choosing r small enough 
to have \Aqq\ < £(n)r ( J" 1 , where Aqq is the second fundamental form of 8Q and e(n) > is a 
dimensional constant to be choosen momentarily. Consider r < ro and p as in the statement. By 
a simple rescaling of the variable by a factor r and a translation, we find B\{p') and new domain 
n' such that \Bi(p')\Q,'\ <h< h and 

dn'nBtip') C {(x,y) € M' 1 " 1 x R : -e(n)|a;| 2 < y < e(n) \x\ 2 } . (4.15) 

Note that, by an analogous computation as above (now /3 := dist(q', r')), we have that 



1) 0Jn-i J [Vl-r 2 - 1 + p - e{n) r 2 j 



r n-2 



dx < h 



< (n - 1) u n -i / [Vl -r 2 -1 + + e(n) r 2 j 



dx. 



One can easily compute (or argue by elementary geometric consideration as before) that h ~ 
^(n+i)/2_ ]\j te that, setting as before v' := —fie n , we have that B\{p' + v') C fl' because of 
(|4.15p . Scaling back to ft, the conclusion hence follows. □ 

4.2. Proof of Theorem I0.lt part I. Here we prove Theorem lO.il (i), (ii) and (hi). 

The proof of (i) follows from the same arguments in Theorem 13. II Indeed, thanks to the local- 
ization in Proposition 14. II we are now in the same position to discard the boundary effects. More 
precisely, by Proposition 12.11 and Lemma 12.21 we know that the E m are uniform A-minimizers. 
Hence, by the regularity in Proposition 12 .51 the sets E m can be parametrized on a optimal isoperi- 
metric ball B ^ 1 by a C 3 ' Q regular function. Hence, we can derive for E m the improved perimeter 
estimate in Proposition 13.21 and use the optimal isoperimetric inequality to conclude (|0.3[) . 

For what concerns (ii), let q E H be a generic harmonic center and let p°P l the center of the 
optimal ball for E m , namely B"^ 1 = B Tm {p'^ t ). Comparing the energy of E m with that of B Tm (q) 
and using that, as shown in the proof of Theorem 13. II it holds |-E m A-B^ p * | < r yr 2 ™ +2 , we get: 

ir%g rm (p°„f) -7& m (?) = F^n(BZ) - F iA B r m {<l)) < ^A B T m ) - F 1>m {E m ) 

\E AB opt \ 2 

< 7 NL(B°£) — j NL(E m ) < 7 ' m n J" J +7C +1 \E m AB°£\ 

<7 2 ^ +3 =7<5ori". (4.16) 

This implies by (|1.18[) in Lemma [TT21 that 

KpT) - Kq) = 9r m (pT) - 9r m (q) + C r 2 n < 5 r» + r 2 m < r 2 m . 

Since the minimum points of h are isolated, from this estimate it follows that p° pt belongs to some 
neighborhood of the harmonic centers. 

Finally, prove (iii) follows readly as in Theorem comparing with the energy of B Tm (p^ 3 *). □ 

5. Stability and exact solutions 



In this section we address the problem of the formation of exact spherical droplets, proving 
assertion (iv) in Theorem lO.il 



IS 
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5.1. Non spherical domains: non existence of critical spherical droplets. In this section 
we show that if Q =^ Br is not itself a ball, the critical points of F ltTn cannot be exactly spherical. 

Proposition 5.1. Let Q C R n be a C 2 bounded open set and assume that it is not a ball. Then, 
Xb t (p) with B rm (p) C fl is not a critical point of F 1 . m . 

Proof. The proof is a simple consequence of a unique continuation argument. Indeed, we show 
that if XB rm ( P ) satisfies the Euler-Lagrange equation (|1.14[) . ()1.15|) . namely 



m 5 



' H 9B rm (p) + 47 V rn = A 

-Av m = XB rm ( P ) -m in tt, 
Vv m ■ v = on dil, 

Ja v m = 0, 

then v rn is a radially symmetric function with respect to p, and hence f2 must be a ball. Assume 
without loss of generality that p = and (|5 . 1 [) holds, and consider the case n > 3 (the two 
dimensional case is analogous). Since Hgs rm = (n — l)/r m , it follows from the first equation 
in (|5.1j) that v m \dB rm = c m G R. Thus, from the uniqueness for the Dirichlet problem for the 
Laplacian, we infer that v m \B rm is radially symmetric and: 

v m (x)= (1 - m) ( J X|2 ~ r " } +c m , for|x|<r m . (5.2) 
An 

Moreover, in D, \ B Tm , v m solves the boundary value problem: 

{Av m = —to in f2 \ B r 

v m = c m , Vv m ■ VQ Brm = (1 ~" )r "' on dB rm . 

Note that also (|5.3[) has a unique solution. Indeed, given vi, t>2 solving (|5 .3[) . w — v\ — v-i solves 



(5.3) 



Aw = in fi \ B r 

w = Vm ■ n = on dB Tm 



(5.4) 



which is extended to a harmonic function in Q setting w = in B Tm , thus implying it; = in 
51 \ _B rm . By a direct computation, the solution of (|5.3[) is given by 

,. M ._ "*(N 8 -d) , . , r2 m C 

m{ h 2n +Cm+ n(n-2) n (n - 2) |x|«- 2 ' 

Therefore, since Vv m ■ 1/ e on 9f2, it follows by the radial symmetry of v m that fHs a ball, which 
contradicts the hypothesis. This gives the desired contradiction and concludes the proof. □ 

Remark 5.2. In particular, in the case of periodic boundary conditions the exact sphere is never 
an equilibrium configuration. 

5.2. Ball domains: uniqueness of a spherical droplet minimizer. In this section we take 
= Br. In this case we show that the ball B Tm is the only minimizer of Fy jm in the regime of 
small mass, thus completing the proof of Theorem 10.11 In order to address this problem, here 
we need to introduce a new ingredient: the stability analysis of the droplet configuration. In 
particular, we will show that the spherical droplet B Tm is strictly stable which will turn to imply 
that it is the unique minimizer of F 1%m . 

Proposition 5.3. Assume £1 = Br C R™, for some R > 0. There exists So > such that, if 
7?"m — ^0 in the case n > 3 and if 7 | log r m | < Sq in the case n — 2, then B Tm is the unique 
minimizer of F 1 , m . 

Proof. The proof of the proposition is made in three steps. 

Step 1. The minimizers E m can be parametrized on B Tm for So small enough. 
To see this, we recall that in the case fi = Br, due to the spherical symmetry, the the origin 
is the only minimum point of the Robin function. Moreover, D 2 h(Q) > Id. To check this, one 
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can either use thethe explict formula for h (see, e.g., [3SJ Chapter IV 5] in the case n = 3, similar 
formulas hold in every dimension): 

or, one can simply notice that R(x, 0) = 2 fl „ , so that D 2 h(0) = £)^_R(0, 0) = WLJ Id fitl . From 
the dchnition of g r in (| 1 . 1 T|) and the radial symmetry of h, it is rcadly verified that g r also has 
minimum in the origin and this minimum is not degenerate as well. From (|4.16[) , we can hence 
conclude that |p^*| 2 <■ ^o r m- Note that, in any dimension n, this implies that 

\p rt\<Sl /2 r m . (5.5) 

This actually leads straightforwardly to the claim. Indeed, for So small enough, there exists 
s < 1 such that, for every point x <G dB rm , B srm (x) n B° pt is a graph over dB rm with small 
Lipschitz constant. Since by (i) of Theorem 10.11 the sets E m are parametrized on dB°^ with a 
graph of small C 1 -norm, this implies that in turns also dE m is a graph on dB Tm . Moreover, the 
C 3 '" regularily is clearly preserved for this new parametrization. 

Step 2. We show now that for So small enough, the ball B Tm is strictly stable. 
Let us recall the second variation for F^ >m . Let E be a stationary point and consider vector 
fields X G C£(fi,R n ) such that 

/ X ■ VEdW 1 = 0. (5.6) 

J BE 

Following [7J Lemma 2.4], for every such field, there exists F : x (— e, e) — > such that: 

(a) F(x, 0) = re for all x G Q, F(x,t) — x for all x € dil and t G (— e, s); 

(b) St := satisfies |£ t | = |£| for every t G (-e,e); 

(c) ^f^k=o = for every x G d£. 

The stability operator for deformations as above is given by (see [3H] in the case n = 2, 3 and 
[1] [16] in the general case) 

F^ m (E)[X] ^Pcr"(E)[X} +1 NL"(E)[X} 

= f (\V 9E (X ■ v E )\ 2 - \A\ 2 (X ■ v E ) 2 ) dW 1 - 1 ^ 

JOE 

+ 87 I [ G{x, y ){X{x)-VE){X(y)- VE )d'H n -\x)d'H n - 1 {y) + 
JdE Joe 

+ 47 f Vvv E (X ■v E ) 2 dU n -\ (5.7) 
Joe 

where \A\ is the lenght of the second fundamental form of dE and v solves (|1.15|) . 

In order to prove the strict stability of B Tm we need to compute (|5.7p on E = B Tm and show 
the existence of a constant Co(n, m, 7) > such that 

<,™(^r m )[X] > Co \\X ■ VB rm \\l> {aBrm y (5.8) 

for every X as in (|5.6j) . We start by noticing that the following inequality holds true: there exists 
a constant ci(n) > such that 

Per"(B r J[X] > ||X • u Brm \\l^ {dBrmV V X| eBrm £ const. (5.9) 

This is a simple consequence of the properties of the spectrum of the Laplace-Beltrami operator 
on the sphere — AQ Brm ■ Indeed, by the decomposition in spherical harmonics (see, e.g., [6l Chapter 
5]), let write X ■ v Br = X)ieN+ fit with fi spherical harmonic in the i th eigenspace of the Laplace- 
Bcltrami operator, i.e. —Ag Br fi = Aj fi and < Ai < A2 < . . . is the discrete spectrum of the 
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operator. Note that i ranges over N + because the harmonic functions on the sphere (i.e., the 
constant functions) are ruled out by the average condition (|5.6[) . If Per" (B Tm ) [X] = 0, i.e. 



dB r 



W 9Brm (X ■ v Br J\ 2 - ^—^ (X ■ u Brm f 



then, by a simple integration by parts and using the orthogonality of the eignespaces, it follows 
that 



As it is well-know, Ai = Hence, Per" \B rm )[X] = if and only if X ■ v Brm = fi is in the first 
eigenspace, that is X\g Brm is constant. Then, the existence of a constant c\ fulfilling (|5.9[) follows 
by the discreteness of the spectrum —A Brm and the obvious scaling property in r m . 

We now recall that the second term in F" m is positive (see [H]); while the last is equal to (for 
E = B rm here v = v m is given in (|5.2[1 ) 

47 / Vv m ■ v Brm (X ■ v Brm f dW^ = - 4 7 (1 = m) rm \\X ■ v Brm f mdB y 

JdB rm n 

Moreover, for a constant vector field X = £ G R", we can compute explicitely NL"(B r7Ti )[C] in the 
following way: 



NL"(i? r J[C] = d2NL &r m (tQ) |^ = d?_ J J { T(\x-y\) + R(x + tC,y + tO)dxdy 



B r „ 

(D 2 R(x,y)((X),(<;X))dxdy>\(f 



t=o 



B r 

^IK-Hli^B^), (5.10) 



because we have already noticed that for f2 = Br the regular part of the Green function has in 
the origin the unique nondegenerate minimum, i.e. Z? 2 i?(0,0) = /ild for some fj, > 0. 

The proof of (|5.8[) is now clear. For a nonconstant X, (|5.8j) follows if 7^ is small enough with 
respect to ci(n). For a constant vector field X (|5.8|) is inferred from (|5.10p . 

Step 3. B r . m is the unique minimizer. The conclusion follows from the fact that the minimizers 
E m are C 2 close to a strictly stable configuration, namely B rm , thus implying that actually E m 
coincide with £>, m . The proof of this fact, well-known for the area functional, can be achieved by 
a carcfull construction of a flow interpolating dE m and dB Tm . Such computations appeared in [TJ 
Theorem 3.9]. In particular, to reduce to this case, let tp m be the paramctrization of dE m /r m on 
dBi, i.e. 

E m = {r m x(l + ip m (x)) : x G dB±}. 

By Proposition ^. 5l and (|5.5j) . for every 77 > we can choose 60 small enough to have H^mllc 3 *" < rj. 
We are, hence, a small perturbation of the fixed stable configuration B\ and [U Theorem 3.9] 
applies. □ 



Remark 5.4. The proof of the previous result becomes trivial if the minimizer E m is such that 
B° E pt is centered at the origin, that is B°g = B Tm . In this case one can indeed drop the quadratic 
term in (|3.5[) and, by the quantitative isopcrimctric inequality, get 

\B rm ^E m \ 2 <jr? n \B rm AE m \ 2 , 

which turns to imply E m = B Tm . 
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Appendix A. On the behavior of the function R in a neighborhood of dil 

In this section we prove the estimates (|1.7[) and (|1.8[) on the regular part R of the Green function 
in (jl.3p . The function R solves 

= M in °> 

^ = Vr x • ^ on c?Sl, 




We introduce the following notation. Since is assumed to have C 2 regular boundary, for x 
in a sufficiently small tubolar neighborhood of <9f2, there exists a unique point xo G 9f2 such that 
dist(x, 9fi) = | a; — xo|. Hence, can hence consider x* G M" \ O such that x* — xq = xq — x and set 
S x := R x + r x ». S x is also characterized by the following boundary value problem: 

AS X = p- in fl, 

\7S X ■ v = (S7T X + S7T X , )-v on 9Q, (A.l) 

J n s x = f n (r x + r x ,). 

The main idea being the estimates are illustrated in the following simple case. Assume that 
G dfl and B 2 H = {x G -B2 : x„ < 0}. Then, by an elementary computation, for every a; G Si, 

J (vr x + vr x » ) • v = on <9fi n s 2 , 
{|(vr x + vr x .)-r/| <i on sn\B 2 , 

Therefore, it follows from (jA.lj) that \S X \ < C. This in turns implies (|1.7|) : namely, there exists 
ro > such that, for r < tq and x,y € fin Bi with r < — x„ < 2 r and |x — y| < r, 

1^)1^1^,(^)1-1^)1. 

Moreover, since for |T X J < |T X | + 1 for every x G B±, (jTTHJ) follows straightforwardly as well. 

The general case of a C 2 bounded domain Q can be deduced by a perturbation of the argument 
above. Let r$ > be such that, for every xo G <9£7, £?2 ro (xo) ^ <^ can be written as the graph of 
a function: namely, up to an affinc change of coordinates, we may assume that xo = and 

B 2ro nil = {(z',t) : t<u(z')}, 

for a given u : B^ C IT" 1 -> K in C 2 (B^) with u(0) = |V«(0)| = 0. In particular, we have 

jn—l 
I ro 

given by 



x = (0, — d) with d = dist(x, dd). Set D := B^^ x [0, 1] and consider the function g : D 



g(z',t) := (vr x ((z',tu(z'))) + VT x ,((z' ,tu(z'))) 



(-tVu(z'),l) 



By definition, g(z', I) = VS^ • i>|gn and g(z',0) = 0. Writing z t — (z',tu(z')), it holds 

d tg (z',t) = u (z') (/-vr x (z t ) + /-vr x ,( Zt )) ■v\ fKh + 
- vr x 0t) + vr x .(zt) . ^ 4= =+ 



^/l + t 2 |Vu(z')l 

<|Vw(z ' )|2 fw. + w..V"k(*). 



l + i 2 |Vu(z')| 2 . 

Since |m(V)| < C|z'| 2 and |Vu(z')l < C|z'|, where C > depends only on | lull (72, one infers that 

\dtg{z',t)\ < (\D 2 r x (z t )\ + \D 2 T x .(z t )\) |z'| 2 + (|vr x | + |vr x 



< 



\x-z t \ n Ix-zA™- 1 ' 



(A2) 
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Note that, for \z'\ < r small enough, 

d + tu(z ] 
d 2 \z'\ 2 



\x - z t \ 2 = \d + tu(z')\ 2 + \zf >t- \ u (z')\ 2 + \zf > | - C \z'f + \A 2 



> 
~ 2 



from which we infer 



\d t g(z',t)\ < , |z,|2 | +, — \ z '\ =; /(/). (A.3) 

' y{ in ~\x-z t \ n |x - 2f |— 1 ~ (d 2 + \z'\ 2 )^ 1 y ' V ' 

It is simple to see that / £ ^(B^ 1 ) for every p£ [1, oo] and 

-. 7T- S ' 1 - 2 d S 

B^ 1 ' ^ {d 2 + S 2 )^ 

f^r t P+n-2 
= d - pin -2)+n-l / « ^ 

{2r if n = 2, 

d _ p( „_2) + n-l r» * P+ '7 2 A < d - P {n-2)+n-l jf „ > 3j (A.4) 
(i+t 2 )^-^ 

Since, for dist(x, dfl) < do, and every z £ dfl fl i?2r- , 

IVS^) • i/|snnB 3ro | = \g{z',l) -g(z',0)\ < / |%?(2:',s)l ds, 

Jo 

we deduce the following bound on the L p norm of VS X ■ v. 

\\VS x -v\\ p LP{aQnB2r) < [ ( f 1 \d t g(z',a)\da) dz> + [ \VS x -v\?dz> 

JannB 2ro \Jo J Jdn\B 2ro 

\d t g{z' , s)\ p ds dz' + C 



anns 2 r Jo 

CO [2r if n = 2, 

< d -p (n-2)+n-l r« f+"- 3 ^ < d - P („-2)+»-l if n > 3. 



< 



Setting /? = (n — l)/p > 0, as a consequence, by L p -regularity theory for (|A.1|) . we get HS^HvyLp 
dodist(x, dVl) 2 ~ n . By the arbitrariness of p and thanks to Sobolev embedding we finally get that 

\S X \ < d^dist(x,dfl) 2 - n . (A.5) 
The proofs of (|1.8j) and (|1.6j) now follows straightforwardly. 

A.l. Proof of (|1.8|) . Fix ro < do as above. Then, if a; belongs to some il ro , then 

\G(x,y)\ < |r(a;,y)| + |i2(a;,»)| 

< \r(x, y )\ + \r(x*, y )\ + \s x \ 

<\T(x,y)\ + l, 

where we used that dist(x, d£l) 2 ~ n < \T(x, y)\ for every y £ H if n > 3, and dist(x, dQ) < 1 in 
n = 2. 

A. 2. Proof of (|1.6|) . Note that, for ro small enough, whenever r < ro/2, r < dist(x, < 2r 
and |y — x| < r, then |r a; *(y)| ~ (j/) | . Then, by (|A.5|) we may assume that do is sufficiently 
small that, for ro < do and x.y as above, it holds 

\R-M\ > \?*-(v)\ - |s*(y)| > |r a (y)| -^dist(.T,a^) 2 -" > |r a ( tf )|, 

and 

\Rx(y)\ < \T x -(y)\ + \S x (y)\< \T x (y)\+d p dist(x,dn) 2 - n < \T x (y)\. 
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